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Abstract. We study stochastic behavior of a single vortex loop appeared in imperfect 
Bose gas. Dynamics of Bose-condensate is supposed to obey Gross-Pitaevskii equa- 
tion with additional noise satisfying fluctuation-dissipation relation. The correspond- 
ing Fokker-Planck equation for probability functional has a solution 7-'({V'(r)}) = 
A/'exp(— {'!/i(r)} /T), where 7if{?/)(r)} is a Ginzburg-Landau free energy. Consider- 
ing a vortex filaments as a topological defects of the fleld i/'(r) we derive a Langevin- 
type equation of motion of the line with correspondingly transformed stirring force. 
The respective Fokker-Planck equation for probability functional 'P({s(^)}) in vor- 
tex loop configuration space is shown to have a solution of the form 7'({s(^)}) = 
J\f exp{—H {s} /T), where A/" is a normalizing factor and H {s} is energy of vortex line 
configurations. In other words a thermal equilibrium of Bose-condensate results in a 
thermal equilibrium of vortex loops appeared in Bose-condensate. Some consequences 
of that fact and possible violations are discussed. 



1 Introduction and scientific background 

Quantized vortices appeared in quantum fluids have been an object of inten- 
sive study for many years (see for review and bibliography the famous book 
by Donnelly Q). The greatest success in investigations of dynamics of quantized 
vortices has been achieved in relatively simple cases such as vortices in rotating 
helium (where they form a vortex array orientated along an axis of rotation) or 
vortex rings. However these simple cases are rather exception than a rule. Due 
to extremely involved dynamics initially straight lines or rings evolve to form 
highly entangled chaotic structure. Thus a necessity of use of statistic methods 
to describe chaotic vortex loop configurations arises. A most tempting way is 
to treat "gas" vortices as a kind of excitation and to use thermodynamic meth- 
ods. One of first examples of that way was an use of the Landau criterium for 
critical velocity where vortex energy and momentum were applied to relation 
having pure thermodynamic sense. More extended examples would be the fa- 
mous Kosterlitz-Thouless description of 2D vortices or its 3D variant intensively 
elaborated currently (for review and bibliography see e.g.[||). 

In the examples above and in many other it is assumed that chaotic vortex 
loop configurations are in a thermal equilibrium and their statistics obeys the 
Gibbs distribution . That supposition is based on fundamental physical principles 
and can be justified in a standard way considering vortex lops as a subsystem 
submerged into a thermostat and exchanging energy with the latter. The role of 
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thermostat in the case of vortices is played by the other excitations (phonons and 
rotons) of an underlying physical field. In the case of Bose- Einstein condensate 
(BEC), which we consider in this paper, that field is an order parameter tp{r,t). 
Vortex lines are just the crossings of surfaces where both real and imaginary part 
of ip{r,t) vanish. Excitations of order parameter (phonons and rotons) interact 
with vortices driving the latter to statistic distribution which in accordance with 
general principles should be the Gibbs distribution. 

It is well known however that the Gibbs distribution can be alternatively 
obtained in the frame of some reduced model like kinetic equations or Fokker- 
Planck equation. That way of course is not of such great generality as a principle 
of maximum entropy, but instead it allows us to clarify the mechanisms of how 
the Gibbs distribution is established and to discuss possible deviations and vi- 
olations of equilibrium state. We choose that way to examine how a thermal 
equilibrium of chaotic vortex loops follows from a thermal equilibrium of BEC. 

2 Langevin equation 

We perform our consideration on the basis of the Gross-Pitaevskii model[^,0|. 
To describe stochastic behavior of BEC let us add to the right-hand side of the 
Gross-Pitaevskii equation the random stirring force ^(x, f). 

Here H {ip}-is the Ginzburg-Landau free energy functional 
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The thermal noise obeys the following fluctuation-dissipation theorem 



(2) 



(C(xi,ii)C(x2,t2)) = ^^1^ <5(X1 - X2) <5(tl - i2) . (3) 

The stochastic problem introduced by has a solution describing thermal 

equilibrium, where probability of some configuration of ^/j-filed is proportional to 
exp{—H Up{r)} /T). This readily follows from the correspondent Fokker-Planck 
equation 0. 

Let's go on to the problem of obtaining of the vortex dynamics appearing 
in BEC, whose own dynamics obeys relations That problem, as well as 

more general problems of that kind has been considered many times by many 
methods]^ starting from pioneering papers by Pitaevskii0 and Fetter We de- 
velop one more method, neither too rigorous nor principally new but convenient 
for purpose of this paper. Details of that method will be published elsewhere, 
here we briefly describe it to an extent to comprehend the put goal. 

As any motion of topological defects the vortex line dynamics is determined 
by the one of the underlying field theory. On the other hand if one ignores a pres- 
ence of other excitations except of the vortices one can say that all motion of 
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BEC is determined by the one of quantum vortices. Thus there is mutual corre- 
spondence and the order parameter iplx, t) can be considered as some functional 
of a whole vortex loop configuration ?/'(x | {s(^,i)}). Temporal dependance of 
field '0 is connected with motion of lines and its rate of change (at some point 
x) is expressed by the following chain rule: 

5V'(x, t) _ f 6^{x I {s(e, t)}) dsie, tl^^,_ 



dt J Ss{e,t) dt 
r 

Vortices in quantum fluids are very slender tubes (except of vicinity of phase- 
transition) and their dynamics is analogous to the one of the strings. That implies 
that we have to aim our efforts at integrating out radial degrees of freedom. It 
can be reached by the following procedure. Let us further multiply equation(|l|) 
by zrjifi ^i^* /Ss{£,o,t), where some chosen point on the curve. Combining the 
result obtained with the complex conjugate and integrating over whole space we 
have 



h f ,3 f A~i di; Si^* ^ A + i di^* 6il^ 
m 



d^x 



+ I dt 5s{£,Q,t) yl2 + l dt Ss{^o,t), 

3 f 6H{^,r) sr 5H{i,,r) h \ 

d ^ -7^^C(x,0 . + -tt-tC (x,t)- 



712 + 1--^ ' '6s{^o,t) ' A^ + 1^ ' ' '<5s(^o,i) 

The first integral in the right-hand side of (^ expresses a chain rule for functional 
derivative dH{s)/Ss{^o,t) where H{s) is the energy of moving BEC expressed 
via vortex line position. Consequently considering them to be very slender tubes 
(which is justified when the radius of the curvature R is much larger of the core 
size ro) and neglecting an energy associated with the core, the quantity H{s) is 
just kinetic energy of the superfiuid flow created by vortices (see e.g. 

PsK f f s'iQs'iO 

""^'^ ^ ^ J J wrm'^'^ ■ 

r r' 

Here s'(^) is tangent vector, double integration is performed along a whole line, 

K is a quantum of circulation equal 2TTh/m. Calculation of functional derivative 
5H{s)/6s{^Q,t) is straightforward and leads to result 

= p,Ks'(eo) X B(eo). (7) 



<5s(eo,t) 



Quantity B(^o) is the velocity of the line element s ( ■fg) expressed by well 
known Biot-Savart law. The terms in the left-hand side of equation (||) (first 
line) can be evaluated in general form by observing that the major contribution 
into integrals appears from vicinity of the vortex filament (see e.g.fO]). Thus to 
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evaluate integral we replace \ s{£_,t)) by ^py{x±) = ^py{s{£_cht) — x) where 
tpv is well studied 2D vortex and integration over d'^x by (P:x.±d£,ci- Functional 
derivative Sijj* / Ss{£^' , t) should be evaluated by a following rule: S^p* /6s{£^' ,t) = 
V_l'0^(x_l)5(^' — ^c/)- Here ^c/ is the label of point of the line closest to point x. 
Using the said above and calculating integrals of squared gradients of tpy (x^ ) 
we conclude that left-hand side of equation (j^) transforms into 

h 2ttPs ■ , ^ , i.^s h 2'KpsaA ■ , , , 

s ( ^o) X s'(^o) + jP-^^iio)- (8) 



Let us now discuss the rest terms of equation (|]) including random force 
C(x, t) (the third line). Consequently considering that the all motion of BEC is 
connected to motion of line, we have to consider Langevin force ^(x, t) as some 
secondary quantity stemming from random displacements of filaments. Con- 
nection between displacements (random) of filaments ds and deviations (ran- 
dom) of (^■0(x, t) may be written in form similar to (|^) with formal substitution 

9-0(x, t)/dt — *■ <5?/'(x, t) and Ss. Taking into account that Sip{x., t) = C(x, t)5t, 
and i5s = ({^,t)6t we conclude that quantities C(x, t) and ({£,,t) are connected 
to each other by a chain rule 

C(x,0=/ ^^^",]g|^j^^" c(e',t)d^ (9) 
r 

That implies that to take into consideration random displacements of line 
we have to change the last term in equation (|^) by the one similar to ( ^ 
with substitution s ( ^q) Ci^Oit). Gathering all terms we obtain a vector 
equation, which can be resolved up to tangential velocity S|| ( ^q) along the 
curve. The latter does not have any physical meaning and can be removed by 
suitable parameterization of the label variable ^. Solving that vector equation 
we arrive at 

« ( ^o) = Y^^^i^o) + ^Y^^^«'(^o) X B(Co) + (m/n)C(Co,i). (10) 



Equation ( |10D describes motion of vortex line in terms of the line itself. It is 
remarkable fact (not obvious in advance that noise C(Coi t) acting on line is also 
additive (does not depend on line variables). 

The last effort we have to do is to ascertain both the statistic properties of 
noise C(?o,i) and its intensity. Shortly, it can be done by comparison of equa- 
tion(||) with substitution s ( ^q) ^ C(Co,i) with the last term of equation (^). 
Clearly the former appeared as result of transformation of the latter. Equating 
them and taking the scalar productions of both parts of the resulting relation 
we arrive at 

(C,,(6,il)C,.(6,t2)) - — ^^^^^i;i^<5(6-6)'^(ti-t2)^,u... • (11) 

PsTT[n/m) 1 -I- yl^cr^ 
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Here and are components of random velocities in rji , 772 directions lying 
in the plain normal to the line. 

Thus starting from dynamics of BEC ( equation (|^)) with the fluctuation- 
dissipation theorem (^) we derive equation ( p^ ) describing motion of vortex line 
in terms of line itself with the additive noise obeying the fluctuation-dissipation 
theorem (pj|). These relations complete a stochastic problem of quantized vor- 
tex dynamics under thermal noise stemming from the one stirring the underlying 
field. In the next section we demonstrate that this problem has an equilibrium so- 
lution given by Gibbs distribution exp{—H {s} /ksT)^ where H {s}- functional of 
energy due to vortex loop (equation ^) and T -temperature of Bose-condensate. 



3 Fokker-Planck equation 

To show it we, first, derive the Fokker-Planck equation corresponding to Langevin 
type dynamics obeyed ( [l0| ) and (|l^). Let us introduce probability distribution 
functional (PDF) 

V{{s{i)],t) = {5{s{0-s{tt))). (12) 

Here 5 is delta functional in space of vortex loop configurations. Averaging is 
fulfilled over ensemble of random force. The Fokker-Planck equation can be 
derived in standard way (see e.g. fl^ ) 



1 + 



B(?) + 

f l)yl(T 



{l + A^)Aa , 



1 + 712^2 



s'(e) X B(C) 



(13) 



Equation ([T3| ) possesses the equilibrium solution in form of the Gibbs distribu- 
tion 7^({s(^)}) = 7Vexp(— if {s} /T), where TV is a normalizing factor. Let us 
show that the first integral term vanishes identically for that solution. To do it 
we exploit relation (^ and parametrization of label variable ^ in which velocity 

s (^0) is normal to the line. Using a tensor notation we rewrite the first term in 
integrand in form (we omit the coefficient and factor exp{—H {s} /ksT)) 



a/37 



SHjs) 



MO 



(5s^(^,t)(5sa(f) 



SH{s) 6H{s) 



The functional derivative ^s^(^)/(5sq(^) oc 5pa therefore all terms vanish due to 
symmetry. Thus the reversible term gives no contribution to flux of probability 
(in the configuration space) equation (|l^), one says it is divergence free. Further- 
more exploiting again relation |^ ) one convinces himself that second (dissipative) 
term in (|3|) and third (due to stirring force) term exactly compensate each other 
(locally) as it should be in the thermal equilibrium. 
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4 Possible violation of thermal equilibrium. 

Thus we have proved that the thermal equihbrium of BEC results in the thermal 
equilibrium of vortex loop. We are now in position to discuss how it can be 
destroyed. Analyzing the proof one can see that the following steps were crucial. 
1. Additive white noise C(x, t) acting on field ip{:x., t) is transformed into additive 
white noise /(f , t) acting on vortex line position s(^,t). 2. Intensity of noise /(^, t) 
expressed by (^l]) is that it locally compensates dissipative flux of probability 
distribution functional in the Fokker-Planck equation ) . 

That observation points out how the thermal equilibrium in space of vortex 
loops can be destroyed for real vortex tangles appeared e.g. in counterflowing 
Hell or formed in quenched superfluids. Being a macroscopical objects vortex 
loops inevitably undergo large scale perturbations generated e.g. by nonuniform 
flow or by action of other vortex loops, randomly placed with respect to the 
studied loop. One more essential source of large scale perturbations might be 
long wave instabilities of vortex filament motion. That type of random action 
drastically differs from (5-correlated in f space thermal noise considered above. 
It obviously cannot compensate dissipative flux of probability, which is pro- 
portional curvature and acts accordingly in small scales. Instead the following 
scenario seems to be realized Due to nonlinear character of the equation of 
motion the large-scale perturbations on an initially smooth filament interact cre- 
ating higher harmonics. They in turn generate harmonics with larger n , where 
K is one-dimensional wave vector arisen in ID Fourier transform (with respect 
to label variable ^) of quantity s(^, t). One can say that an additional curvature 
created by large scale stirring force propagates in region of small scales. In real 
space that corresponds to entangling of vortex loop and creation of vortex line 
segments with large local curvature. Then dissipative processes come into play, 
their role now is that they remove from the system harmonics with very large 
wave vectors k. So the Kolmogorov cascade-like solution with a flux of a curva- 
ture in space of k is established and, as a result, the stochastic distribution is 
far from equilibrium. That scenario is quite similar to the one which is realized 
in classical turbulence or in the so called weak (wave) turbulence. A difference 
is that if in case of the wave and classical turbulence) an exchange an energy 
between harmonics is realized, whereas in our case there is exchange a curvature 
(k''s(k)s(— k)). 

The scenario described above is especially relevant in the so called low- 
temperature superfluid turbulence case, when the normal component is very 
small and the usual (considered here) dissipation monotonically increasing with 
the curvature is absent. Some recent experiments and numerical simulations 
show that vortex tangle in Hell decays at extremely low temperature (about 1 
mK), where dissipation due to normal component is negligibly small. Obviously 
some other strong mechanisms of dissipation must take place. It can be e.g. 
emission of phonons and rotons from speedy moving parts of line or just col- 
lapse of "hairpin" segments of the filament. That mechanisms are concentrated 
on very small scales (of order of the vortex core size) or equivalently in region 
of very large wave numbers k. Therefore regions of the pumping and the sink 
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of an additional curvature are greatly remote in k space and a system must be 
essentially noncquilibrium. 

One more reason of a violation of thermal equilibrium might be a reconnec- 
tion of lines. In real vortex tangle consisting of many loops the vortex filaments 
undergo frequent collisions and rcconncctions. Just after reconnection there ap- 
pear kinks on the curves disappearing later on. From mathematical point of 
view a kink on the curve can be described as a discontinuity of tangent vector 
s'(^,i) which has Fourier transform of type s'{k) oc k~^. Thus the reconnec- 
tion processes supply the selected curve with discontinuities having a spectrum 
s{k) oc k"^ . Taking into account a random nature of vortex line collision the 
reconnection processes can be modelled, in some measure as a random stirring 
of filament with spectrum of type {({k)(^{—k)) oc k~^. Remember now that the 
establishing of a thermal equilibrium requires that the random force correlation 
function is proportional cx 5(^1—^2) and, consequently the spectrum docs not 
depend on wave number k, {({k)(^{—k)) = const . Therefore the colored noise 
{C{k)C{—k)) oc k~'^ coming from reconnection processes can also lead to the 
Kolmogorov type noncquilibrium state. To clarify which of mechanisms forming 
noncquilibrium state prevails in real vortex tangle one has to investigate much 
more involved problem. 

This work was partly supported under grant N 99-02-16942 from Russian 
Foundation of Basic research. 
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